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Abstract 

We introduce NeedATool (Needlet Analysis Tool), a software for data analysis based on needlets, a 
wavelet rendition which is powerful for the analysis of fields defined on a sphere. Needlets have been 
applied successfully to the treatment of astrophysical and cosmological observations, and in particular 
to the analysis of cosmic microwave background (CMB) data. 

Usually, such analyses are performed in real space as well as in its dual domain, the harmonic one. 
Both spaces have advantages and disadvantages: for example, in pixel space it is easier to deal with 
partial sky coverage and experimental noise; in harmonic domain, beam treatment and comparison 
with theoretical predictions are more effective. During the last decade, however, wavelets have emerged 
as a useful tool for CMB data analysis, since they allow to combine most of the advantages of the two 
spaces, one of the main reasons being their sharp localisation. 

In this paper, we outline the analytical properties of needlets and discuss the main features of the 
numerical code, which should be a valuable addition to the CMB analyst's toolbox. 
Subject headings: methods: data analysis, numerical, statistical, cosmology: observations, cosmic 
background radiation 



Over the last two decades, the detailed analysis of 
cosmic microwave background radiation anisotropics has 
been fundamental in determining the global properties 
of our Universe and its evolutions. The cosmological 
concordance model encodes into few parameters the va- 
riety of processes we observe in the local universe as 
well as those occurring at very large scales. Such pa- 
rameters have been measured very prec i sely by sev- 
eral CMB experiments (Mather et al. 1992; Sm oot et al. 
119921 : Ide Bernardis et al.l 12000b iKomatsu et al.l l2010f ). 
and such measurements will be further refined with the 
next generation of cosmological experiments. CMB data 
analysis is very demanding - both in terms of compu- 
tational power required and sophistication of the neces- 
sary techniques - due to the complexity of the datasets 
and to the high degree of accuracy one wants to achieve. 
An overview on the techniques recently applied to as- 
trophy sical data analysis can be found in lPesenson et all 
([201(1 This holds in particular for cutting-edge experi- 
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ments such as the ongoing Planc k satellite 1 . 

Oy e r the last decade, wavelets dFreeden and Schneider l 
1998; Antoine and Vanderghevnst 1999; McEwen et al. 
I2QQ6L 12007b ISanz et al.l 12006b IStarck et all I2QQ6L I2009D 
have emerged as a very powerful tool for CMB 
data analysis; a very incomplete list of references 
should include testing for non-Gaussianity (|Vielva et al.l 
120041: ICabella et all 120041 ). foreground subtraction 
Han sen et al.l 120061 ), point sourc e detection (ISanz et al.l 



componen t separation ([Moudden et al.l 120051: 
Starck et al.ll2006j ), polarisation analysis ([Cabella et al.l 



20071 ) and many others. The reason for such a strong 



interest is easily understood. As it is well-known, CMB 
models are best analysed in the frequency domain, where 
the behaviour at different multipoles can be investi- 
gated separately; on the other hand, partial sky cov- 
erage and other missing observations make the evalua- 
tion of exact spherical harmonic transforms troublesome. 
The combination of these two features makes the time- 
frequency ^ocjdisj^ip^ wavelets most valu- 
able. See Starck and Bobin (2009) for a recent review on 

1 http://www.rssd.esa.int/SA/PLANCK/docs/Bluebook-ESA- 
SCI(2005)l.pdf 
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this topic 2 and IWiaux et al.l ([20081 ) and for a numerical 
implementation. 

Recently, a novel approach to spherical wavelets 
has been introdu ced in the statistical literature by 
iBaldi et al.l (|2006[ ), adaptin g tools proposed in the func- 
tional analysis literature bv lNarcowich et ah! (|2006[); the 
first a pplication to CMB data is due to iPietrobon et al.l 
(|2006[ ), where needlets are used to estimate (cross- 
) angular power spectra in order to search for dark energy 
imprints on the correlation between larg e scale structures 
and t he CMB (jSachs and Wolfe] \l96% . IGuilloux et all 
(|2007| ) investigate the effect of differ ent window f u nctions 
in needlets constructions; whereas IBaldi et all (|2009a| ) 
provide further mathematical results on their behaviour 
for partially observed sky-maps. Needlets have been 
applied to angular power spectrum estimation in the 
prese nce of noise (|Fav et al.l 120081 : iFav and Guillouxl 
120081 ). as well as to the estim ation of the bispec- 
trum (|Lan and Marinuccil l2008f ); the latter tool has 
been applied to the WM AP 5-year da t a rel ease by 
IPietrobon et~aT1 (|2009D and lRudiord et al.l (|2009D to con- 
strain the primordial non-Gaussianity parameter. The 
bispectrum formalism has been further e xploit ed by 
IPietrobon et~aTl (|201Q[ ) and iRudiord et all (|2010f ). who 
addressed the sky asymmetry issue analysing respectively 
the properties of 3-point correlation func tion, and the 
primo rdial non-Gaussianity parameter. iCabella et al.l 
(|2010f ) developed the bispectrum estimator including 
the marginalisation over the possible foregroun d resid - 
uals in the CMB m aps, while iDelabrouille et al.l (|2009f ): 
iGhosh et al.l (j2010[ ) produced a foreground component 
separation algorit hm. The an a lysis of directional data 
are described in IBaldi et all (f2009bf ). Finally, the 
needlet formalism has bee n extended to the the polarisa - 
tion field, as discu s sed by IGeller and M arinucci (2008); 



Geller and Ma^el (120081. l2009h: IGeller et al.l (j20O 
Geller et al.l (j2QQ9h : IGeller and Marinuccil (|2010l ). 



The aim of this paper is to describe a numerical code, 
called NeedATool (Needlet Analysis Tool), which has al- 
ready been used in several of the above mentioned analy- 
ses. We first provide a discussion of the needlet formalism 
in Sec. [I] and [2j We then describe a viable implementa- 
tion of the code in Sec. [3] and we add our concluding 
remarks in Sec. 2J 

1. NEEDLETS FRAME 

Needlets enjoy several features which are not in gen- 
eral granted by other spherical wavelets construction. 
He re, we recall some of t hese properties an d refer 
to iMarinucci et al.l (j2008f ): lLan and Marinuccil (j20Q8h 
for a comprehensive mathematical discussion. Com- 
plementary math e matica l analyses can be f o und i n 
Geller and Mavelil (l2QQ7h: lLan and Marinuccil (|2QQ9h : 
Mavelil (|2008l ): IGeller and Mavelil (|2008l . I2009D . In par- 
ticular, needlets 

a) do not re ly on any tange nt plane approximation 
(compare iSanz et aLll2006[ ). and take advantage of 
the manifold structure of the sphere; 

b) being defined in harmonic space, they are 
computationally very convenient, and natively 

2 A mutliresolution package for data analysis and compression 
is also available at this url: http://jstarck.free.fr/mresol.htm 



adapted to standar d packages such as HEALPix 1 
(jGorski et al.ll2005h : 

c) they allow for a simple reconstruction formula (see 
Eq. [5]), where the same needlets functions appear 
bot h in the direct and the inve rse transform (see 
also iKerkvacharian et al.l (|2007| )): 

d) they are quasi-exponentially (i.e. faster than any 
polynomial) concentrated in pixel space, see Eq. [6] 
below; 

e) they are exactly localised on a finite number of 
multipoles; the width of this support is explicitly 
known, controlling the power encoded in each mul- 
tipole range (see Eq. [I]); 

f) needlets coefficients can be shown to be asymptoti- 
cally uncorrelated (and hence, in the Gaussian case, 
independent) at any fixed angular distance, when 
the frequency increases. 

We first recall that the spherical needlet function is 
defined as 




) £ Y im (j)Y £m (Z jk )- (1) 



where 7 and £jk are directions on the sphere, Y^ m is a 
spherical harmonic function, with Y ' i m identifying its 
complex conjugate, and b(x) is a filter function defined 
for x G [1/B,B], which the entire needlet construction 
relies on. Here, we use {Qk} to denote a set of cubature 
points on the sphere, corresponding to frequency j; and 
Xjk denotes the cubature weights. In Fig. [T] the needlet 
profile as function of the angle between 7 and £jk is shown 
for the choice B=2, j=8. 

Nccdlct shape; B = 2, j = 8 



Fig. 1. — Needlets in pixel space. B = 2, j 
of the angle between 7 and £jk 



as a function 



Intuitively, needlets should be viewed as a convolu- 
tion of the projection operator Y?m=-£ Y z™{7)Y£m(€jk) 
with a suitably chosen window function b(-). The needlet 
frame construction strongly relies on the spherical har- 
monic decomposition which represents the mathemat- 
ical environment for the derivation of the fundam en- 
tal properties of needlets (|Lan and Marinucci l2009f ): in 

1 http://healpix.jpl.nasa.gov 



particular, the existe nce of a reconstruction f o rmula , 
as first pointed out b y Freede rTand Windheuserl (|1997h : 
IFreeden and Schneiderl (11998J). M ore details can be 
found in iFreeden and Maier (2002) . A similar approach 
has been followed bv lStarck et al.l (2006. 2009), who im- 
plemented wavelets, ridgelets and curvelets built directly 
on the sphere, both for scalar and spin-2 fields. All 
these sets of functions are examples of frames on the 
sphere: they are over-complete and redundant but admit 
a well defined backward transformation. For this reason, 
frames are not only suitable for a multi-frequency anal- 
ysis of a signal, as any other wavelet implementation, 
but also for data compression, denoisin g algorithm and 
component separation, as descr ibed bv iMoudden et all 
([20051 ): IStarck and Bobinl (|2009h . 

Besides spherical harmonic decomposition, needlet 
properties strongly depend on the the filter function fr(-), 
which controls the angular scale span covered by each 
needlet and ensures that needlets enjoy quasi-exponential 
localisation properties in pixel space. Formally, we must 
ensure that: 

i) b 2 (-) has support in and hence b(^j) has 
support in £ £ [B^.B^ 1 ] 

ii) the function &(•) is infinitely different iable in 
(0,oo). 

iii) we have 

oo p 

^6 2 ( — )^lforalH>5. (2) 

3 = 1 

It is immediate to see that property (i) ensures the 
needlets have bounded support in the harmonic domain; 
property (ii) is the cruci al element in the derivat ion of the 
localisation properties ([Narcowich et al.l 12006): finally, 
property (iii) is necessary to establish the reconstruc- 
tion formula (Eq. |5j). Functions such as b 2 (-) are called 
partitions of unity. 

There are of course many possible constructions sat- 
isfying (i-ih); indeed an interesting theme of research is 
the derivation of opt imal windows satisfy ing these three 
conditions (compare iGuilloux et al.ll2007[ h although the 
choice of &(■) is expe cted to exert second-order effects on 
the final estimates (|Lan and Marinuc ci 2009). An ex- 
plicit recipe for the construction of &(■) is given in Sec.[3l 
Very recently an extended study on how needlet prop- 
ert ies depend on the filte r functions has been conducted 
by iScodeller et al.l (|201Q[ ) . Interestingly, the authors ex- 
plore a peculiar construction able to mimic the Spherical 
Mexican Hat Wavelets based o n mathematical study of 
IGeller and Ma^el (12001 l2009h . 

Needlets coefficients are hence given by 

Pjk = V^jk^Z 6 i^jj Yl a trnYi m (£ jk ). (3) 

l ^ ' m=-£ 

In Fig. [2] we show the needlet coefficients of WMAP 5- 
year temperature map for the specific choice B = 2 and 
j = 4. A remarkable aspect of this construction is that 
the needlet coefficients can be represented easily as a 
Mollw eide projection in t he HEALPix pixelization frame- 
work (jGorski et al.ll2QQ5h . the most widely used tool for 
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visualization and analysis of CMB maps. This makes 
dealing with needlets particularly handy, since it is easy 
to implement a needlet analysis code which exploits pre- 
existing HEALPix routines. 

QVW Needlets Coefficients j=4 




Fig. 2. — Needlet coefficients of the combined Q, V, W map 
at the resolution j = 4. The B parameter i s fixed to 2. No- 
tice as the anomalous bright spots found bv lPietrobon et aTl 
(|2008T ) are clearly visible. 

It is very important to stress that, although the 
needlets do not make up an orthonormal basis for square 
integrable functions on the sphere, they do represent a 
tight frame. In general, a tight frame on the sphere is a 
countable set of functions {ej} such that, for all square 
integrable functions on the sphere / G L 2 (S 2 ), we have 

so that the norm is preserved. Here (f,g) means the 
scalar product, or more properly the projection between 
the functions / and g. Of course, this norm-preserving 
property is shared by all orthonormal systems; how- 
ever, frames do not in general make up a basis, as 
they admit redundant elements. They can be viewed 
as the closest system to a basis, for a given r e dund ancy, 
seelHerna ndez an d Weiss! (|1996l ). IBaldi et al.l (|2006l ) and 
iBaldi et al.l (|2009a[ ) for further definitions and discussion. 
In our framework, the norm-preserving property trans- 
lates into 

j,k 1=1 

where 

Ci = 2ZTT^' a£m ' 2 

m 

is the angular power spectrum of the map T(ff). Identity 
3] has indeed been verified by means of numerical simu- 
lations and implicitly provides the correct normalization 
for needlets, Ajk- It is basically a consequence of the pe- 
culiar partition-of-unity property of &(•) (Eq. [2]). Eq. [4] 
is related to a much more fundamental result, i.e. the 
reconstruction formula 

W) = ^2Mjk(i), ( 5 ) 

3,k 

which in turn is a non-trivial consequence of the ana- 
lytical properties of the bej functions. As mentioned 
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before, the simple reconstruction formula of Eq.[5]is typ- 
ical of tight frames but does not hold in general for other 
wavelets systems. 

The following quas i-exponential loc a lisatio n property 
of needlets is due to iNarcowich et al.l (j2QQ6h and moti- 
vates their name. For any M = 1,2,... there exists a 
positive constant cm such that for any point x G S 2 we 
have 

1^(7)1 < „ | D , ; w m ■ (6) 



(1 + Bi arccos((x, Cjk)) 



We recall that arccos((x, is just the natural distance 
on the unit sphere between the points (x, £jfc). The mean- 
ing of Eq. [6] is then clear: for any fixed angular distance, 
the value of ipjk (7) goes to zero quasi-exponentially in the 
parameter B. This clearly establish an excellent localiza- 
tion behaviour in pixel space. Note that the constants 
cm do depend on the form of the weight function &(•), 
and in particular on the value of the bandwidth param- 
eter B; typically a better localisation in multipole space 
(i.e., a value of B very close to unity) will entail a larger 
value of cm, that is, less concentration in pixel space 
for any fixed j. The resulting trade-off in the behaviour 
over the harmonic and real spaces is expected: smaller 
values of B correspond to a tighter localisation in har- 
monic space (less multipoles entering into any needlet), 
whereas larger values ensu re a faster decay in real space 
(jLan and M arinucci 2009j). 

In iBaldi et al.l ( 2006[ ), another relevant property of 
needlets coefficients was discussed, namely their asymp- 
totic uncorrelation at any fixed angular distance, for 
growing frequencies j. More explicitly, at high frequency 
needlets coefficients can be approximated as a sample 
of identically distributed and independent (under Gaus- 
sianity) coefficients. Also, in view of Eq. [3j for full sky 
maps and in the absence of any mask we should expect 
the theoretical correlation to be identically zero whenever 
|.7*i — .72 1 > 2. This has be en indeed numerically verified 
bv lMarinucci et al.l ([20081 ). 

The probabilistic properti es of the needlet co efficients 
f3jk have been established in lBaldi et a l. (2009a); in that 
paper, it is shown that for any two (sequence of) pixels 
€jk,€jk' such that their angular distance is larger than a 
positive e, for all j, we have 



cm 



M-l 



for all M = 1,2,3, 



(7) 



thus proving wavelets coefficients are asymptotically mi- 
correlated as j —> 00 for any fixed angular distance. Eq.[71 
can then be seen as the statistical counterpart of Eq. [6l 
These properties are the basis for the large success of 
needlets as a CMB toolbox, in particular when dealing 
with masked datasets. 

2. NEEDLETS ESTIMATORS: 2- AND 3- POINT 
CORRELATION FUNCTIONS 

Having introduced the spherical needlet frame, and 
recalled the main properties which make needlets per- 
form extremely well in a wide number of applications 
to 2-dimension fields on the sphere, we now briefly de- 
scribe some important statistical techniques largely used 
in CMB data analysis. 



2.1. (Cross-) Power Spectrum 

After computing the needlets coefficients f3jk from a 
2-dimension map (e. g. the CMB or source count map), 
we can use Eq. [31 to build a (cross-)correlation estimator 
in wavelet space, as: 



1 si «j 



(8) 



where N p i x (j) is the number of pixels (e. g. in the 
HEALPix scheme 7Vpi x = 127V s 2 ide ) with I and J denoting 
the two different maps. The theoretical prediction for f3j 
can be computed from the expected Cj J as: 



^ J = E 



(21+1) 

47T 



t 



(9) 



where we recall Cj J = (<4 m a^) is the (cross-) angular 
power spectrum. 

Pj provides then an unbiased estimator for the (cross-) 
angular power spectrum within the needlets framework. 
The analytic relation between f3j and Ct underlines few 
more advantages in using needlets. Indeed, it makes 
extremely easy and straightforward dealing with beam 
profiles and experimental window functions, which have 
to be taken into accou nt when analysing real data (see 
iPietrobon et al.l ([20061 )). The duality which needlets em- 
bed, namely the localisation both in pixel and harmonic 
domain, allows also to ch a racter ise t he noise proper- 
ties (see iDelabrouille et all ([2009) and IPietrobon et al.l 
(|2009L [2010D for direct applications to WMAP 5-year 
data.) 

Computing the 4-point correlation function, it can be 
easily shown that the analytical expression for the disper- 
sion of the estimated cross-correlation power spectrum in 
needlet space is: 



A/3" 



\ 



£ 



(21+1) 
16tt 2 



which, of course, must be only taken as an approxima- 
tion when dealing with real data, where window func- 
tions, noise and partial sky coverage have to be taken 
into proper account. 
It is important to stress that Eq. [8] generalises into 



pix 



J 

32 k 1 



(10) 



which describes th e needlets coeffic i ents c ovariance and 
it has been used in IPietrobon et al.l ([20 08) to determine 
the degree of anomaly of a few hot and cold spots found 
in the CMB temperature map. 

We have shown that the needlets formalism may be 
suitable for the problem of angular power spectrum es- 
timation from a CMB map (and therefore, indirectly, to 
the estimation of cosmological parameters). In particular 
the application of needlets to the WMAP 3-year data led 
to interesting constrain ts on the dynamics of dark energy 
(jPietrobon et al.ll20Q6h and to the measure of the differ- 
ence in power between the two estimates of the power 
spectrum computed on the north and south CMB skies 
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([Pietrobon et al.ll2008[ ). A detailed discussion on the ap- 
plication of needlets to pow er spectrum estimation can 
be found in iFav et all (|2C 



2.2. Needlets Bispectrum 

In the previous section we described how needlets can 
naturally be applied to the estimation of the 2-point cor- 
relation function and how, thanks to the reconstruction 
formula (Eqs. [2] and 2]), it relates to the usual angular 
power spectrum. It is easy to extend the formalism to 
the higher order correlation functions. 

Here, we focus on the 3-point correlation function, 
which plays a crucial role in CMB data analysis to de- 
tect any departure of primordial fluctuations from the 
Gaussian statistics, a smoking gun for non-standard in- 
flationary models. We next briefly review the properties 
of the needlet bispectrum and how it relates to the spher- 
ical harmo nic bispectrum. An exten si ve discussion is 
provid e d in lLan and Marinuccil (120081): IPietrobon et al.l 
(120091): [Rudiord et al.l (|2009D : IPietrobon et al.l (|201Q[ ): 
iRudiord et al.l (|2010D . 

The needlet bispectrum is defined as follows: 



$ 313233 — jy < Pj 1 fc flj 2 kflj 3 k 



(11) 



{31)1X32)1(32,) 



(24 + 1) (2*2 + 1) (24 + 1) 



47T 







Bp 



where 



is the estimated bispectrum, averaged over ra^s. Sj 1 j 2 j 3 
can be seen as a binned bispectrum, a smooth and 
combined component of the angular bispectrum. The 
bispectrum is supposed to be vanishing for a Gaus- 
sian d istribution. Standard inflation mechanism (iGuth 



1981 



iSatol 119811 : iLindd 119821 : lAlbrecht and Steinhardt 
1982|) predicts a tiny non-Gaussianity in the cosmo- 
logical perturbations: this is why a great effort has 
been spent to measure a bispectrum amplitude differ- 
ent from zero in the CMB data, which would provide 



an extraordinary h a ndle on the early universe physics 




ter ms of the non-linear parame te r /nl (see for exam - 
ple iKomatsu and Spergel (1200 lh : iBartolo et all (|2004f ): 
iSmith and Zaldarriagal ([2006[ )). 

One of the key properties of needlets, is that the sum 
of the squared filter functions in harmonic space, 6^-, is 
1 (see Eq. [2j). This means that even if we group multi- 
poles and each needlet peaks at a certain multipole range 
the total power is conserved. Therefore, the needlets 
power spectrum analysis can be in principle done with 
any choice of the parameter B, being the total power 
conserved and only the correlation and localisation prop- 
erties affected by different width of bij. We incidentally 



notice that any wavelet function defined as the difference 
between the square roots of the scaling functions at two 
different resolutions satisfies Eq. [2j but the uncorrelation 
properties are in general not granted, since they are de- 
termined i ^y_the_sJm2eof thefilter function fe^(-). See for 
example iScodeller et al.l (|2010f ) for a detailed discussion 
on needlet filters. 

This does not hold any more when the the cubic power 
of the filter functions contributes to the estimator used, 
which is indeed the case of the skewness expression and 
more generally of the bispectrum one. This fact is dis- 
played in Fig. [3] where we plot the sum of square and 
cube of the filters functions. The not uniform sampling 
of the multipoles for a n-power estimator suggests that 
the choice of the B parameter is crucial for the analy- 
sis and must be driven by the insight on the range of 
multipoles to be probed. 



b, sum 



1.2 r~ 
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0.8 -\ 
0.6 - 
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Fig. 3. — Solid line sum of the bj; dot-dashed line sum of 
the b\. While the former is equal to 1 for the entire range of 
multipole, the latter is not. 



3. NEEDATOOL: A NUMERICAL IMPLEMENTATION OF 
NEEDLETS 

In the previous sections we have discussed the main 
needlet properties, which are indeed strictly related to 
the analytical properties of the filter function in har- 
monic space, b(£/Bi). We now give a specific recipe 
for the construction of &(•), and describe the main fea- 
tures of "NeedATool" Needlet Analysis Tool - a numer- 
ical code which computes the filter functions and the 
needlet coefficients according to the algorithm which 
has been applied in the several cosmological analyses 
(IPietrobon et al.l l200l 120081 120091 1201(1 ICabella et all 
I20ia . 

We give a step-by-step pr ocedure, as impl e mente d in 
NeedATool and described in iMarinucci et al.l (|2008l ). 



1. Construct the function 

exp(-T372) , 



/(*) 



f ex] 



-1 < t < 1 

otherwise 



It is immediate to check that the function /(■) 
is C°° and compactly supported in the interval 

(-1.1); 

2. Construct the function 



ip{u) = 



J\f(t)dt 
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The function tjj(-) is again C°°; it is moreover non- 
decreasing and normalised so that ip(— 1) = , 
V(l) = l; 

3. Construct the function 

( 1 if < t < ± 

<p(t) = { ^(i-^r(^-i)) if i< t <i 

[0 if t > 1 

Here we are simply implementing a change of vari- 
able so that the resulting function (/?(•) is constant 
on (0,£> _1 ) and monotonically decreasing to zero 
in the interval (£? _1 ,1). Indeed it can be checked 
that 

!_Jfl_ (4 _I )= r 1 for t=\ 
B-V B> \ -1 for t = l 

and 

V(^)=^(l) = l 
¥>(1)=^(-1) = 0; 

4. Construct 

fc 2 (0 = v(|)-v(0 

The expression for b 2 (-) is meant to ensure that 
the function satisfies the partition-of-unity prop- 
erty of Eq. O Needless to say, for &(£) = 
f p 1 1/2 

<(p(|) - | we take the positive root. 

Incidently, we notice that property |4] is crucial in al- 
lowing for the reconstruction formula, and it is shared, 
although within a different setup, by the implement a- 
tion described by iStarck et al.l ([20061 l2QQ9h who makes 
use of the HEALPix software package too. In Fig. [4] we 
show the set of filter functions in £ space for the choice 
B = 2. They result in a homogeneous binning in log£, 
whose power sum to 1, a crucial property which needlets 
properties rely on. 



b- filter 

1.2 I — 1 




10 100 1000 

1 

Fig. 4. — Filter function in ^-space which the needlet con- 
struction relies on. Set computed for B = 2. 



NeedATool 2 computes the needlet filter functions in 
harmonic space and the needlets coefficients, given a 



http: / /www. fisica.uniroma2.it/~pietrobon/dp_webpage_eng. html 



set of parameters. The cod e is based on the public 
available HEALPix package 3 (Go rski et al.ll2QQ5h . which 
has to be downloaded and installed separately (version 
2.10 or higher required). The software is composed of 
two programs "synneed" and "ananeed"; following the 
HEALPix structure the former deconvolves a given map 
to obtain a needlet frame, the latter reconstructs the orig- 
inal map, if a needlets basis is given as input. Both the 
programs, as well as the routines necessary for NeedA- 
Tool will be fully integrated in the next HEALPix re- 
lease. Both the programs accept the same parsing file in 
which the fundamental parameters are provided by the 
user. A list of such parameters is given in Tab. [TJ The 
maximum number of multipoles (4nax) and the B pa- 
rameter are required. Those given, the codes computes 
the maximum j necessary to keep all the information in 
the map. The N s ^ e of the needlets coefficients is then 
determined according to the relation £ max < 2_V S id e - The 
filter functions bgj are computed by default, while it is 
possible to choose whether computing the needlet coeffi- 
cients, which actually is the most time-consuming part, 
by setting the keyword " compute_needlets" . It is neces- 
sary then to specify the map and its resolution. A sky 
mask can be applied filling the "maskfile" variable. The 
last two keywords set the output files. The input map 
of ananeed has to be the needlet output file created by 
synneed. 

The cubature points Xj = {£,jk} k= i 2 are assumed 
to coincide with the pixelization of the unit sphere S 2 
provided by HEALPix , with A^de such that / max = 
< 2 A^ide (with [•] denoting the integer part 
and B > 1). The cubature weights, Ajk are given by 
1/Apix, with Npi X given by 12 • N^ de . We then com- 
puted the /3jk coefficients for each k position given by 
the HEALPix scheme evaluating the projection opera- 
tor, namely the product of YgmX eii'm' for each pair 
of pixels ^jfc, by means of the HEALPix software pack- 
age. The code is very fast, and it can be run on a laptop. 
For a low resolution map side — 256), it takes a few 
seconds, while it scales according the HEALPix scaling 
laws for higher resolutions. 

We present an example of the reconstruction power of 
the needlet frame. We produced a CMB map consistent 
with t he WMAP 5-year best fit cosmological power spec- 
trum (|Komatsu et al.l f2QQ9h up to ^ max = 500 by using 
the HEALPix toolbox. 

We processed this map through the needlet pipeline 
extracting first the needlet coefficients by applying syn- 
need, then reconstructing the map using ananeed. We 
repeated this procedure both for the whole sky case and 
in presence of symmetric sky cuts, 15°, 5°, 1° and 0.1°. 
The results are shown in Fig. [5] (left column) together 
with the percentage error due to the procedure (right 
column). The error due to the forward and backward 
transformation is smaller then 0.01% except for few pix- 
els in the whole sky case. The result worsens reaching 
1% when a broad mask is applied, basically because the 
needlet coefficients at very large scale are affected by the 
presence of the mask. This effect is actually expected, 
since needlets work well at small angular scales, namely 
when the filter function b(£/B^) groups several multi- 

3 http://healpix.jpl.nasa.gov/ 
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TABLE 1 

Parameters required by the codes "synneed" and "ananeed" . 



Parameter 


synneed Ananeed 


healpix_dir 


/usr/local/Healpix_2.13a 




500 


B 


2.0 


compute_needlets 


T 


mapfile 


input/lcdm_map_lmax500.fits needlet_2.00_Nj009.fits 


mapnside 


256 256 


maskfile 


input / sky _cut_l_256_ring. fits input / sky _cut_l_256_ring. fits 


bl2_root 


bl2 


need_root 


needlet recmap 



poles. The mask effect may be reduced by either fine 
tuning the B parameter, or deconvolving the mask effect . 
Reconstruction error is indeed wavelet frame dependent: 
a better performance can be achieved when the wavelet 
is defined as the difference between two scaling functions 
at two different resolutions instead of the difference be- 
tween the square root of the s caling functions. Th is is the 
implementation discussed by iStarck et al.l (|2QQ6[ ) which 
leads to an exact reconstruction. 

To give another estimate of the error due to the re- 
construction, we computed the angular power spectrum 
of the original map, C^, as well those of the reprocessed 
ones, Of, and compared them. The ratio is displayed in 
Fig. [6] (left panel) for the four analysed sky cuts applied. 
The full sky reconstruction shows an excellent agreement, 
the difference between the power spectra being of the or- 
der of 10 -4 ; when a mask is present, the reconstruction 
causes an error of few percent at very low angular scales, 
which decreases at small scales. 

The non perfect reconstruction in presence of miss- 
ing information is expected. Therefore, it is interesting 
to compare the needlets performance to the one of the 
spherical harmonics implemented by HEALPix . We ex- 
tracted the a£ m from the same CMB realisation for the 
five different sky coverages considered by using anafast; 
then we produced a CMB map out of the pseudo-aim 
using synfast. We finally analysed the resulting map 
computing its power spectrum. The difference with re- 
spect the power spectrum computed through needlets is 
shown in the right panel of Fig. [6j The spectra agree 
very well regardless the applied mask. This confirms the 
high performance of needlets in reconstructing a field in 
the sphere. As a further figure of merit, we show in 
Fig. Q the percentage error between the two maps recon- 
structed using spherical harmonics and needlets for the 
broadest mask (15°). The agreement is striking: a large 
scale pattern appears at the level of < 1%. 

We further investigate the mask effect, focusing on how 
the presence of a sky cut couples the spherical harmon- 
ics coefficients. As it is well known, in the case of partial 
sky coverage the spherical harmonics do not form an or- 
thonormal basis. This can be formalized as follows: 



T(j)Y > m / d^ = ^ atmKtt' 

ira 



(12) 



(13) 



(14) 



/3jk coefficients as: 

I'm' 

rR (7) = E E K Jj) a<w Y, b {w) Yemmu ' mml - 



j t'm' 



where O is the observed region and TCwmm' is the cou- 
pling matrix. In the needlet framework, this affects the 



(16) 

As for the spherical harmonics decomposition, the map 
reconstructed through the needlet pipeline is computed 
from a£ m ; in addition, the mask shows indirectly its effect 
on the filter functions b(£/B^). 

To visualize this effect, we plot in Fig. [8] the coupling 
matrix Kwmm' fc> r a full sky case (upper panel), in which 
it actually becomes a Kronecker's 8 function, fe/ mm /; for 
a 15° symmetric sky cut (middle panel) and for the con- 
volution ^2jb(£' /Bi)K,££'mm'W which appears in 
Eq. [161 Since the numerical evaluation of the coupling 
matrix up to £ = 500 is a severe computational challenge 
and we expect the higher multipoles to be marginally 
affected by a partial sky coverage, we computed Eq. HI] 
on a smaller range of spherical harmonics in the interval 
£ £ [0, 20] on the HEALPix pixelization at N side = 16. 
The presence of the mask translates into off-diagonal 
terms in Fig. [8j each multipole is coupled to its second 
neighbor, whereas the coupling with the first neighbor 
is vanishing because of parity. The coupling is indeed 
not negligible up to the fourth second neighbor. The 
lower panel in Fig. [8] clearly depicts the properties of 
needlets and explain the error we obtain in the recon- 
struction. The correlation matrix results a superposition 
of blocks, each of them corresponding to a given j: each 
block partially overlaps only the first neighbor, as a con- 
sequence of the compact support of the b(£/B^) func- 
tions. (This is shown by the white regions we observe 
in the plot). Moreover, the off-diagonal elements are less 
powerful compared to the auto-correlation terms and the 
correlation extends to a lower number of multipoles. At 
very large scale however, a lack of power is present in 
the diagonal terms which affects the global reconstructed 
map. This can be seen in the right panel of Fig. [5j where 
the error in the reconstructed map follows a peculiar pat- 
tern which results symmetric as a consequence of the sky 
cut. 

The additional coupling due to the mask present in the 
needlet construction arises from the definition of needlet 
itself, Eq.[U since the needlet is built from the projection 
operator. However, needlets provide a natural way to 
counterbalance this effect. At large scale, the needlets 
are not sharply localised, and the coefficients turn of to 
be non- vanishing also in the masked region. Taking into 
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Fig. 5. — Left: reconstructed maps for the several sky cuts analysed (from top to bottom full-sky, 15°, 5° and 1°); right: percentage 
error on the reconstructed map. 
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Fig. 6. — Left panel: absolute value of the difference between the power spectra extracted from the original map and the reconstructed 
one for several symmetric sky cuts: full sky (red solid line), 15° (yellow solid line), 5° (light blue solid line), 1° (green solid line) and 0.1° 
(blue solid line). The black solid line shows the CMB angular power spectrum for the original map. Right panel: comparison between 
the reconstruction using synfast-anafast combination and the needlets pipeline for the same sky cuts. The input angular power spectrum 
contains information up to t = 500, but we computed the power spectrum of the realization map up to t = 700 to check the leakeage from 
low multipoles to high ones: this explains the edge at i = 500 in the plotted power spectra. 

4. CONCLUSIONS 

In this paper we have introduced NeedATool, a public 
software for the analysis of datasets on the sphere based 
on the needlet framework. The software is particularly 
useful for the analysis of CMB data, as shown by its 
successful application to, e.g., the WMAP dataset. The 
needlet construction differs from other wavelet renditions 
due to the distinctive properties of the filter functions bgj , 
which translate into a sharp localisation in pixel space 
and excellent properties of non-correlation among the 
functions of the set. This aspect is crucial when building 
estimators for CMB data analysis as we discussed exten- 
sively: therefore, needlets are a very promising tool for 
high-accuracy cosmological experiments. 

Fig. 7. — Percentage error between the two maps reconstructed 
using spherical harmonics and needlets. A symmetric sky cut of 
15° is applied. The agreement is impressive. 




account this leakage, and including in the reconstruction 
procedure the coefficients inside the mask, it is possible 
to reduce the effect of the coupling between needlets. 
This is the actual procedure implemented in the code, 
which leads to the results discussed in this section. The 
residual effect due to the presence of the mask shows up 
in the large scale pattern highlighted in Fig. [71 
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